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The three-band model with the 0-0 dh'ect hopping near to unit fihing is considered. We present 
the general procedure of reduction of this model to the low-energy limit. At unit filling the three- 
band model in the charge-transfer limit is reduced to the Heisenberg model and we calculate the 
superexchange constant. For the case of the small electron doping the three-band model is reduced 
to the t — J model and we calculate electron hopping parameters at the nearest and next neighbors. 
We derive the structure of corrections to the t — J model and calculate their magnitude. The values 
of the hopping parameters for electron- and hole-doping differ approximately at 40 %. 
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1. Introduction 

Since the discovery of electron-doped su- 
perconductors there is growing interest in 
structure and properties of these systems 0, ^ 
Q . The electron- and hole-doped high-tempera- 
ture superconductors (such as Nd2-a:Cea:Cu04 
and La2_a;Sra,Cu04) both have the Cu02 planes 
with the same structure constant. In spite of 
the similar structure the experiments show the 
strong asymmetry in the properties of these 
compounds. Thus, the critical temperature for 
Nd2_i.Cea;Cu04 is only 22 K that in contrast 
with T, = 40 K for La2_^Sr^Cu04 0. Also, 
the doping concentration which destroy the AF 
order for electron- and hole-doped systems are 
~ 15% and 2 ~ 3% respectively ||, |. To 
study of such asymmetrical phenomena might 
help to illuminate the nature of superconduc- 
tivity in these compounds. 

Since the structure of the Cu02 plane is the 
same for these compounds, one can expect that 
well known three-band Hamiltonian should 
be suitable for both of them. 

In this paper we investigate one electron 
over unit filling in a Cu02 plane. We start from 
the Emery Hamiltonian |^ with the 0-0 direct 
hopping. On the basis of our general approach 
to the low-energy reduction 0, we obtain the 
effective single-band Hamiltonian. The equiv- 
alence of this single-band Hamiltonian and the 
t — J model is discussed and the second-order 
corrections to the t — t' — J model are derived. 
The effective parameters of the t — t' — J model 
for the electron are calculated in the realistic 
region of Emery model's parameters. Effective 
hopping is less than in the case of hole-doping. 
Similar conclusion was obtained in the work 
by Zhang and Benneman 0], but they used 
the perturbation theory over the parameters 
tpd/Ud, tpd/{ep — ed) which does not work as it 
was argued in our previous works 0, |^. 

In Sec. II we represent the three-band model 



in more suitable terms for reduction to the 
single-band model. We use a three-step pro- 
cedure. At the first step we introduce the sym- 
metrical and antisymmetrical oxygen operators 
At the second step we separate out the one- 
site Hamiltonian and get its solution. At the 
third step we represent the primary Hamilto- 
nian in terms of the Hubbard operators which 
reflect the structure of the solution of the one- 
site problem. In Sec. HI we briefly consider 
the three-band model at unit filling. Then, we 
apply the Schrieffer- Wolff transformation [llO 



for getting the J-term of the low-energy Hamil- 
tonian. In Sec. IV we get the effective electron 
hopping parameters and correction to them. 
Also, we discuss the accuracy of the t — t' — J 
models for the electron-doped system. Section 
V presents our conclusions. In Appendix some 
details of our consideration are given. 



2. Three-band Hamiltonian 

The three-band model with the direct 0-0 
hopping is given by the following Hamiltonian 

i 



La 



+U,J2dl^d,^dl^di^ + H' , (1) 
I 

where (dia) creates (annihilates) the hole 
in a Cu d^-^-y-^ state at site /, {pmp) cre- 
ates (annihilates) the hole in a O Px{y) state at 
site m, ed and are the energies of Cu and 
O levels respectively, Ud is a intrasite Coulomb 
repulsion at the copper site. 

The hybridization term H' is given by 

H' =t Yl (dlPma + H.C.) 
<lm>a 

-tp (PmaPm'a + H.C.) , (2) 

<mm'>a 
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where <lm > denote the nearest-neighbor Cu H' = 2i(: ^ ^U'{dta1i'a + H.c.) (5) 

and O sites, < mm' > denote the nearest- <ii'>a 

neighbor O sites. In Eq. (|) we follow to the -tp ^ {fJ'ii'i'liaQi'a - oL^Ua) 

sign convention of Ref.]!!], II^]. <"'>" 

As it was firstly noted by Zhang and Rice ~^^u'{Qi 

+ H.C)}, 

T3|, it is conveniently to use the orthonormal- , r r^i m ■ . ^ 

— \ . . r 1 where the exphcit form of the coeiricients A^i', 

ized oxygen states on the tour oxygens around , i , ■ t r -r^ /m\ /f^ 

. . ^ \m u /U//', I///' can be obtained from Eqs. (yj), (g), ^) 

a Cu site. In our previous work |7| we have 

represented the Hamiltonian (P, (||) in terms {A, fi, v}uf = {A, fi, z/}(l - 1') 

of orthogonal symmetrical and antisymmetri- = ^{A n i^}kexp(— zk(I — I')) (6) 

cal oxygen cluster states k 



Qi 



E(l + 7k)"'/'(cos(A;,/2)pkx ^^^h 



+ cos{ky/2)pi^y) exp(-ikr. 



Ak =(l + 7k)'/' 



i/2/_„„„^i, /oN„, (o^ /ik = 8cos2(A;^/2)cos2(A;j,/2)(l + 7k) 



qi =E(l+7k)"'/'(-cos(V2)Pk. (3) 

k i^k = 4 cos(fc2:/2) cos(/cy/2)(cos^(/cj;/2) 

+ cos(A;j2)pky) exp(-2kr;) , _ cos^(A;y/2)(l + 7k)"^/^ . 

where summation in Eq. (D is produced over These coefficients decrease rapidly with increas- 

the Brillouin zone with 7k=(l/2)(cos(k^a)+ ing |1 - r|. The values of A, jj and u for small 

cos(^a)) and pkx,y are the Fourier image of |i _ i'| are given in Table I. One can easily get 

Pmx,y that the uqq = and, therefore, the mixing 

E/ N /.N of o and q at the same sites is absent. The 

expf— ikr^,) . (4) r i tt -i • /ft\ • • i 

y transformed Hamiltonian (g) is equivalent to 

the three-band Hamiltonian (|l|), (^). 
The physical reason for introduction of qi and ^j^j^g ^^e Hamiltonian (D into local 
qi states |jl3| is as follows. The hole at the ^^^^ hopping parts 
copper can hope only at the symmetrical com- 
binations of the oxygen states. If we separate -f^Zoc = ^ d X! ^Ux^ia + (cp — /^o^p) X! ^ta^ia 
out the oxygen states which interact strongly 

with the copper states, we can solve the prob- +(ep + /xo^p) X! ^ta^ia. 
lem of determination of the low-energy two- 

hole states or vacuum states at the background +Ud X] d^di^^dl^dii (7) 
of the one-hole states (spins). The states qi, qi 

(^) are independent at different sites and are -|-2tAo ^(fi^q'ia + -ff.c.) , 
ortonormalized. They are very convenient for 

solving this problem. -f^/iop = 2t ^ ^ii'idux^i'a + H.c.) 

The original Hamiltonian ([1|), in terms 

of qi, qi takes the form —tp^if'ii' iliali'a — QiaQi'a) 



ll'a 



Ho =edJ2 (^tadia + epJ2{qt^qia + qtjia) +'^ii'{qtqi'a + H.c)} , 

l,a l,a 

j^jj^y2 (i+ dii hereafter sum over Z, I' denotes I ^ I'. One can 

I '"^ see that hybridization term in Hi^c (0) includes 
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only the symmetric oxygen state in agreement 
with Zhang and Rice |T^ . The direct 0-0 



hopping, once taken into account, does not mix 
local states with the opposite symmetry. 

For the case of unit filling there is one hole 
per unit cell. Therefore, one can introduce the 
set of one-hole cluster state 

\do, >= d+\0 > , \q^ g+|0 > , 
>=g+|0>, (8) 



and rewrite Hioc in this terms 



m 

la 



La 



qq 

la 



+2t\,{Xf^ + H.c.)} 



(9) 



where 



ab 
la 



Wla X b, 



la\ 



(10) 



where Ri = {A^ + At'^XlY/\ A = {A-fiotp)/2, 
and e/,g = -Ai t Ru Ai = (A + fiotp)/2, 
A = ep — ed- As discussed in Ref. 0, we as- 
sume that at unit filling the groundstate of the 
Cu02 plane is the low \fa>- states on each 
cluster (the hole on the copper with admix- 
ture of symmetrical combination of the hole on 
the nearest four oxygens) with a virtual tran- 
sitions at neighbors. Such transitions give the 
superexchange interaction in the second order 
of perturbation theory. 

The hopping part of the Hamiltonian (j^) 
contains transitions between vacuum, one- and 
two-hole states at the different sites. For a de- 
tailed consideration of this subject see Ref. [0]. 



3. Superexchange interaction 

Here we consider only the terms of -ff^op 
that are relevant to intersite interaction 



is the Hubbard operator at the site I, a = 
±^ is the spin projection. It is convenient 
to introduce the Hubbard operators because 
they form the natural basis for description of 
one-site states. If we introduce also the non- 
diagonal Hubbard operators we can simply ex- 
press all operators in our Hamiltonian (|^ in 
their terms. Diagonalization of if/^^ (^ is per- 
formed for each site independently. After di- 
agonalization Hl^^ is given by 



Hi 



hop - 2^ i'w,l3 \-^l 
ll'al3,y 

+xrxr-} , 



0,//3 



(13) 



^loc 



H {^f^fa + (^gXla 



where y is the set of two-hole states, which 
are five singlets and three triplets (Ref. 0), 
F^,'^|^ are matrix elements of the Hamiltonian 
for transition from {/;, ///} to {yi, 0//} states. 

By applying the Schrieffer- Wolff transfor- 
mation to (|l^) (see Appendix and [0, |10|) one 
can get 



with 



L. 



+ {ep + f^otp)X!!} 



\fa >= U\dc, > -V\qa > 
\ga >= V\da > +U\qo, > 



Hj = Y,{Ju'SiS,+Ya>NiNv] 



(14) 



(11) 



;i/2)a,^Xf^ 



(12) 



U 
V 



{{Ri 



A)/2Rif^ 
A)/2i?i)^/2 



H 1 — 

the expressions for the Jui and Yui constant are 
given in Appendix. 

Thus, we established that at unit filling the 
groundstate of the Cu02 plane is the system 
of local spins which interact antiferromagneti- 
cally. Since Jw decreases rapidly with increas- 
ing |1 — I'l, we hold only the nearest-neighbor 
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term in Eq. ([T^ ) and receive the Heisenberg terms of the t — t' — J Hamiltonian. But the 
Hamiltonian. One can check that y<//'> — expression for the hole-singlet hopping is more 
— i J<///> and so complicated 



<ii'> ^ 
J = J<w>. (15) 

The second term in the Hamiltonian (|15[) can 
be essential when the system goes away from 
unit filling. 

4. Electron-doped system 

The doping of the Cu02 plane by the elec- 
tron is equivalent to removal of a hole from one 
of the clusters. Hence, in the hole picture the 
electron is 'zero' or vacuum state. If the hole 
(spin) is taken away from cluster, the neighbor- 
ing holes (spins) can hope to the empty site. 
Thus, the mechanism of the movement of the 
electron is different from the movement of the 
added hole. The electron moves as a 'hole-in- 
hole', whereas the hole moves as a local singlet. 



Tj'i, = 2tXu, ( V2Ui U-WiV)x 

x{V2ViV -WiU) (18) 
-tpl2a'{V2V,V-W,U)y2 , 

where U and V were determined in Eq. (|T^), 
Ui, Vi, Wi are the components of eighenfunc- 
tion of the local singlet. 

Thus, we mapped the three-band model for 
one electron over unit filling to the t — t' — J 
model. The Hamiltonian (^),(|I6]) in terms of 
the Hubbard operators can be rewritten in a 
more usual form 

Ht-j =t ^k^Cl'a 

<ll'>,a 

+J J2 i^i^i' - 7^/^/') . (19) 
<//'> ^ 

Si = (l/2)c/"o-Q , fii = {cl ■ Cl) . 

Here c^,Qq, are the electron creation and an- 
nihilation operators. 



4-1. Hopping Hamiltonian 

From expression for H^op (0) with defini- 
tions ( p!^ and (ID) one can easily get the zero- 
order electron hopping Hamiltonian 



IVa 



with 



(16) 



(17) 



The electron hopping Hamiltonian (|IB]) coin- 
cides with the hole hopping Hamiltonian de- 
rived in 

works ||. They correspond to the t — t' — ... 



4-2. Second-order corrections 

With the help of Scrieffer- Wolff transfor- 
mation we shall obtain the second-order cor- 
rections to the t — t' — J Hamiltonian. The 
corrections to the first hopping parameter will 
be small and the t — J model may be valid. The 
correction to the other hopping parameters are 
not relatively small and, therefore, models with 
the hopping at the neighbors farther than near- 
est must be more complicated than the simple 
t — t' — . . . type. 

First of all we shall obtain correction to the 
energy of the electron. Such correction arises 
due to the virtual process of hopping of the 
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electron at the neighbors with transition in ex- 
cited state and back. Corresponding Hamil- 
tonian is given in Appendix. The resulting 
second-order correction is 



5He = -Y.MwX'^'^Nv 



(20) 



with 



Mw 



^11' J 



\rpg 12 

+tpfiii>UV. 



(21) 



The Hamiltonian 6He ( pOD coincides com- 
pletely with 6He in Ref. 0, but the expres- 
sion ( piD for the Mill is sufficiently more simple 
than for the hole one. 

Corrections to the direct hopping parame- 
ters can be obtained by transformation of the 



Hamiltonian (|13D. The result is 



6Ht 



'l2lnl'al3,y Tinli{XyXi 
'^^yX'i'^^ Xl""'^ (cr a/jSn) } 



(22) 



where y is the set of two- hole states, Xy = 1/2 
for singlets and Xy = 3/4 for triplets, 2;^ = 1 
for singlets and Zy = —1/2 for triplets. The 
expressions for T^^^, are given in Appendix. As 
it was shown , the corrections to the t — t' — 
J Hamiltonian have a complicated nature and 
depend on the filling (iV-term) and spin state 
of the neighbors (S-term). The corrections for 
the hole single-band model and the electron 
one are quite different. 



4.3. Quantitative results 

We have argued that our general ap- 
proach to the low-energy reduction of a three- 
band model is essentially more correct than 
the other one, because we construct the per- 
turbation theory over the ratio of the effective 



hopping parameters between the different local 
states to the energy gap between them. This 
ratio for the case of CUO2 plane is of the or- 
der of 0.1 . All values of the effective single- 
band Hamiltonian have been derived for ar- 
bitrary parameters of the three-band model, 
without any assumption about a relationship 
among t, A and Ud ■ In our calculation we 



take Ud - 
Up = Vpd 



8.2 eV, t = 1.4 eV, 



0.7 eV, 



according to the different papers 
. Also, we have taken experimental 
126 meV IITTI and have deter- 



value of J ^ ^ 

mined selfconsistently charge-transfer gap to 
be A = 5.07 eV. In Table II we present the fol- 
lowing effective parameters of one-band model: 
the hopping parameters to the first, second and 
third neighbors; the contribution of the direct 
0-0 hopping; corrections to hopping and en- 
ergy ; value of the ratio ti/ J. 

We obtain for the electron hopping t*^ = 
ti ~ 0.39 eV that in 1.36 times less than the 
hole hopping |0] t'* ~ 0.53 eV . The role of the 
direct 0-0 hopping for the electron-doped sys- 
tem is lower than for hole one. Actually, the 
admixture of the 0-state for one hole on a clus- 
ter is small, unlike for the hole singlet, where 
added hole mainly on the oxygen. Thus, the O- 
O hopping contribution to the hopping at the 
nearest-neighbors for electron (tpi/tiY ~ 20% 
and for hole (tpi/ti)'' ^ 36%. 

As for the validity oft — J,t — t' — J and 
other models to the electron-doped system, our 
conclusions are similar to the case of the hole- 
doped system. Thus, if we hold only the hop- 
ping at the nearest-neighbors, corrections to 
the t — J model on a ferromagnetic background 
is near to 0.13% for electrons and 2.5% for 
holes, and t — J model is valid. If we are in- 
terested in the next-neighbors hopping, it is 
necessary to include the correction (P^ ) to the 
effective Hamiltonian. This correction has the 
complicated structure and does not reduce to 
the simple direct hopping. 
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5. Conclusion 

We have studied tlie electron doping of the 
Cu02 plane in the framework of three-band 
model. By applying our general procedure of 
the low-energy reduction to the electron-doped 
system, we conclude that the t — J model is 
valid as well as to the hole-doped system. We 
established some asymmetrical properties. Ef- 
fective hopping parameter for the electron is 
in 1.36 times less than the hole one. Correc- 
tions to the single-band electron and hole t — J 
models have the different magnitudes. Such 
corrections are important for transitions at the 
next-nearest neighbors. 

It is doubtful that the difference in the pa- 
rameters of one-band model may lead to the 
drastic difference in behavior of hole-doped and 
electron-doped systems with doping. The 
strong asymmetrical dependence of the critical 
concentration for disappearance of antiferro- 
magnetizm is probably connected with the na- 
ture of antiferromagnetic ordering in this sys- 
tems In the hole-doped systems antiferro- 



magnetizm has a quasi-two-dimensional char- 
acter and is associated with Cu spins. In the 
electron-doped systems the contribution of Nd 
or Pr spins in the antiferromagnetic ordering 
is essential and three-dimensional effects are 
more important. 
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Appendix: The second-order corrections 
to the hopping process 

For obtaining such corrections we use the 
Schrieffer- Wolff transformation 

H^H = exp(-5)ifexp(5), 

S+ = -S. (23) 

The first order generator of transformation and 
the second-order correction are 

[Ho,S] = -H\ 6H' = {1/2)[H',S], (24) 

for this model Hq = Hioc (Eq.(0)). For the 
derivation of J-term H' = Hhop (Eq.(|l^)). Ex- 
pressions for the Jw and Yiu integrals are given 
by 

\pV |2 



E 



X7, 



ey + eo - 2e/ 

\fw\ 



^ ^ + eo - 2e/ 



(25) 



where y is the set of two-hole states at a cluster, 
ey are there eighenenergies, eo = is the energy 



of vacuum state. Values =1/2 for 



singlets and Xy = —1, Zy = 3/4 for triplets. 
Matrix elements F^^, are obtained in Ref. [0]. 
One can check that the major contribution to 
J and Y is originated from the matrix element 
of the transition to the lowest singlet and so 



<U'> — i^<iv>- 
The terms of Hhop (0), which lead to the 
energy correction, are given by 



ll'a 

. (26) 



Corresponding generator of the transformation 
is 



E 



rp9 



ll'a ^9 



(27) 
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where |/q, > and \ga > are the one-hole states References 



(|l^) with their energies ej and eg. Thus, the 
second-order correction to the energy has the 
form (|20|). 

Corrections to the direct hopping (^21) can 
be obtained from the transformation of the 



Hamiltonian 

by 



-'-InV 



Magnitudes Tf^i, are given 



+ eo-2e/ 



(28) 



The expressions of Ff^ were obtained analyti- 
cally and computed for definite values of pa- 
rameters of the three-band model in Ref. 0. 
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Table I. 

The values of the coefficients A(I — I'), — I'), — 1') as functions of (1 — 1') = rax + my 



,ni ^n,m ^m,n l^n,m l^m,n ^n,m ^m,n 



, 


0.9581 


1.4567 


0.0 


1,0 


0.1401 


0.5497 


0.2678 


1 , 1 


-0.0235 


0.2483 


0.0 


2,0 


-0.0137 


-0.1245 


0.0812 


2,1 


0.0069 


-0.0322 


0.0609 


2,2 


0.0035 


0.0231 


0.0 



Table II. 

The first three hopping parameters, the contributions of the direct 0-0 hopping to them, second-order 
corrections to them on a ferromagnetic background, second-order corrections to the energy, and the 
ratio ti/J. 



neighbor direct hopping direct 0-0 hopping Corrections 
number, n tt/f, (eV) t;fj , (eV) St^ff, {%) 



1 -0.3896 -0.0765 0.125 

2 0.0180 -0.0345 122.5 

3 0.0480 0.0173 24.2 



Correction to Ratio 
the energy (cV) -0.0931 \ti/J\ 3.081 
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